Strong correlation, electron-phonon interaction and critical fluctuations: isotope 
eff"ect, pseudogap formation, and phase diagram of the cuprates 
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Within the Hubbard-Holstein model with long-range Coulomb forces, we revisit the charge- 
ordering scenario for the superconducting cuprates and account for the presence or the absence 
of a formed stripe phase in different classes of cuprates. We also evaluate the mean-field and the 
fluctuation-corrected critical lines for charge ordering and we relate them with the various pseudogap 
crossover lines occurring in the cuprates and we discuss a mechanism for their peculiar isotopic 
dependence. Considering the dynamical nature of the charge-ordering transition, we explain the 
spread of T* and of its isotopic shift, obtained with experimental probes with different characteristic 
time scales. 

PACS: 71.10.-w, 74.72.-h, 71.45.Lr, 71.28.+d 



I. INTRODUCTION 



The phase diagram of the cuprates, hke, e.g., 
La2-a;Sra;Cu04 (LSCO), is characterized by three dis- 
tinct regions which are reached by varying the doping 
x: (i) the overdoped region, where the system displays a 
markedly metallic character, which can be reasonably de- 
scribed within the standard Fermi-liquid picture; (ii) the 
region around optimal doping, where the superconduct- 
ing critical temperature Tc is maximum, and no other 
energy scale besides the temperature T is seen in various 
experiments, so that the metallic phase above Tc strongly 
deviates from the Fermi-liquid behavior; (iii) the under- 
doped region where, on the contrary, there is an abun- 
dance of energy scales, and an even stronger violation of 
the Fermi-liquid behavior. These are related to a doping- 
dependent temperature T*^, below which a suppression of 
the density of states at the Fermi level is revealed in static 
measurements and to a doping-dependent temperature 
T* below which various experiments reveal a pseudogap 
Ap. T* strongly depends on the characteristic time scale 
of the probe used to measure the pseudogap properties. 
Other energy scales of the underdoped phase are the su- 
perconducting gap A and the critical temperature T^. 

The generic partition of the phase diagram into three 
qualitatively quite distinct regions with markedly differ- 
ent behaviors and the absence of any energy scale but 
the temperature in the intermediate (optimally doped) 
region, suggest the presence of a quantum critical point 
(QCP) near optimal doping ||l|-^. If this is the case, 
the identification of the optimally doped region with the 
region where quantum critical fluctuations are present 
naturally accounts for the peculiar non-Fermi-liquid be- 
havior. Various proposals have been made to specify the 
broken symmetry phase to be associated with the under- 
doped region, ranging from a circulating-current phase 
, to a change in symmetry of the superconducting order 
parameter The scenario that we have been propos- 
ing along the years ||,^ is related to the occurrence of an 



instability for charge ordering (CO), along a line Tco{x) 
ending in the QCP near optimal doping [^-12|. This 
proposal of a (dynamical) CO instability finds support 
from neutron scattering experiments in Ne-doped LSCO 
JlSf and YBa2Cu307_x systems Q, where charge do- 
main walls (stripes) have been detected. 

In this paper, by extending previous results, we ac- 
count for the various pseudogap-formation temperatures 
and give a reason why formed stripes may be present in 
some classes of cuprates and be absent in other classes. 
We also address the issue of the anomalous isotopic ef- 
fects of Tc and T* in the cuprates and provide a novel 
mechanism based on the relevance of CO fluctuations. 



II. THE CHARGE-ORDERING SCENARIO AND 
THE CUPRATE PHASE DIAGRAM 



The above CO scenario found one of its possible real- 
izations within a minimal model of strongly correlated 
electrons interacting with phonons, i.e., the one-band 
Hubbard-Holstein (HH) model in the presence of long- 
range Coulomb forces ||l0| , [ri]] . In this model there is a 
hopping term t representing the electron kinetic energy, 
a (large) on-site Hubbard repulsion U, representing the 
strong electron-electron (e-e) correlation, a coupling g 
between a local lattice distortion and the local electron 
density, a phonon frequency ujo, characterizing the rel- 
evant lattice distortion, and a long-range e-e Coulomb 
repulsion Vc- The simplification with respect to models 
taking into account the complex structure of the cuprates 
was adopted for two main reasons: first of all it makes the 
model manageable despite the formal complications re- 
lated to the treatment of strong e-e correlations via the 
slave-boson large- expansion formalism. Moreover it 
allows to concentrate on the relevant aspects providing 
both a non-Fermi-liquid behavior and a strong pairing 
mechanism. Although the relevant features related to 
quantum criticality should be generic, different cuprates 
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should of course show diflFerent aspects related to their 
specific structure. 

Two major results were obtained within the HH model 
First of aU, at mean field, one obtains the T = 
phase diagram reported in Fig. 1 with a line of the 
electron-phonon (e-ph) coupling g as a function of the 
doping X. This line marks a second-order transition for 
the onset of CO characterized by an order parameter 
Pqc representing the microscopic charge modulation at 
a wave vector psf . Each point on the curve g{x) is 
a QCP, which is the T = end point of a critical line 
T^q{x). Secondly, charge collective fluctuations mediate 
a doping- and temperature-dependent singular scattering 
among the quasiparticles 



r(q,c^)^c/ 



V 



(1) 



where 17 is a residual repulsion, V is the strength of the 
singular interaction, 7 is a characteristic time scale 

of the charge fluctuations, and ^q"^ is the inverse square 
correlation length, which measures the distance from crit- 
icality. This interaction diverges ai — and at a critical 
wave vector q^ when vanishes. This occurs at T = 
when the doping is reduced down to a mean-field critical 
value and at finite T, when approaching the mean- 
field critical line T^q{x). As we discuss in Sec. II B, the 
CO fluctuations beyond mean field shift the mean-field 
critical line Tqq{x) to a lower value Tcoix), ending in 
the fluctuation-corrected QCP at a;c (< a;"), and modify 
accordingly the expression for the coherence length ^. 

The singular part of the interaction ([|) (with Co ^ C if 
fluctuations are taken into account) is attractive both 
in the particle-particle and in the particle-hole chan- 
nels. Therefore, the quasiparticles feel an increasingly 
strong attraction by approaching the CO critical line. In 
the particle-hole channel, this interaction can produce a 
pseudogap due to the incipient CO. At the same time in 
the particle-particle channel, the strong attraction can 
lead to pair formation even in the absence of phase co- 
herence JTgI . Therefore the pseudogap- formation temper- 
ature T*{x) closely tracks the underlying CO transition 
line. 

We point out that the assumption of an infinitely large 
Hubbard repulsion U and the large- expansion do not 
allow for a correct description of the spin degrees of free- 
dom jT^. However, at large but finite U the above sce- 
nario is not expected to drastically change, except that 
the modulation of the charge profile within the CO state 
should favor antiferromagnetic ordering in the charge- 
poor regions. This state with charge and (enslaved) 
spin modulation mimics the stripe phase observed in the 
cuprates, when approaching it from the large-doping side, 
where charge degrees of freedom play a major role. En- 
slaved nearly critical antiferromagnetic fluctuations with 
a characteristic wave vector q^^ ~ (tt, tt) are expected to 



mediate an interaction similar to Eq. (P in the particle- 
particle channel, which becomes repulsive in the Cooper 
channel. This additional contribution is relevant to prop- 
erly reproduce the main features of the single-particle 
spectra ||l^,|l9|. The d-wave symmetry of the supercon- 
ducting order parameter, which is already favored by the 
residual repulsive term in Eq. (|l|) is further stabi- 
lized by the antiferromagnetic fluctuations. 
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FIG. 1. Phase diagram e-ph coupling g vs doping x of the 
single-band infinite-(7 Hubbard-Holstein model with nearest- 
neighbor hopping t — 0.5 eV, next-to-nearest-neighbor hop- 
ping t' — —{l/6)t, phonon frequency ujq — 40 meV, and 
in the presence of long-range Coulomb forces with strength 
Vc ~ 0.55 eV. See Ref. for a detailed description of the 
model. An empty circle on the g axis marks the minimum 
value 3min allowing for CO. 



A. The zero-temperature phase diagram 

The separation between the homogeneous phase and 
the CO phase at T = 0, reported in Fig. 1, is marked 
by a line having a minimum value gmin, below which no 
CO (or stripes) could be formed as a static phase. Of 
course strong quantum fluctuations mediating the singu- 
lar scattering (^, which leads to pairing and non- Fermi- 
liquid behavior in the normal phase, would be present 
both above and below the quantum critical line. 

The crucial point in obtaining the CO instability is 
that strong electronic correlations reduce the effect of the 
homogenizing electron kinetic-energy term, thus favoring 
the possibility of an instability induced by a residual in- 
teraction, which in our model is provided by the e-ph in- 
teraction. The relevant role of the lattice in driving the 
instability is, indeed, strongly suggested by the peculiar 
isotopic effect of T* and ^ and by recent EXAFS 
experiments pH , which identified interesting features as- 
sociated to the lattice structure. Specifically a micros- 
train £ of the Cu-0 bonds was measured in terms of the 
deviation of the Cu-0 distance with respect to a refer- 
ence distance do, which introduces a mismatch in the 
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lattice between the Cu02 layers and the rock-salt layers. 
Actually there is a close resemblance between the exper- 
imental {e,x) and the theoretical {g,x) phase diagrams 
indicating that the Cu-0 microstrain and the e-ph cou- 
pling g of the simplified HH model are strictly related. 
From the microscopic point of view it is indeed quite 
natural that a lattice contraction in the Cu02 planes can 
enhance the effective coupling between the electrons and 
the ions. Therefore the comparison between experiments 
and theory allows to draw the conclusion that a one-to- 
one monotonic relation g = g{s) is likely to exist. 

Once this general framework is settled, one can move 
to identify more detailed possible scenarios. In particular 
the order of the homogeneous-metal to charge-ordered- 
metal transition is of obvious relevance. The simplest 
possibility is that in the real materials this transition is 
of the second order. In this case, at T = 0, the different 
microstrains determined by the rock-salt layers directly 
correspond to different e-ph couplings and are reflected in 
the similar phase diagrams. Of course this simple theo- 
retical picture is not necessarily realized in the quite com- 
plex real materials where i) anharmonic effects and/or 
ii) additional non-ordering fields can partially or entirely 
transform the second-order transition into a first-order 
one. In this last case, the above scenario maintains its 
validity, as long as this transition is weakly first-order. 
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FIG. 2. The phase diagram of the cuprates according to the 
CO-QCP scenario for LSCO (a) and Bi2212 (b). The solid 
line is the m-f critical line Tqq{x) ending at T = in the m-f 
QCP at x^- The lowest dashed line in panel (a) marks the 3D 
critical line Tco{x) in the presence of fluctuations, ending in 
the QCP at Xc- The dot-dashed line in panels (a,b) indicates 
the "dynamical instability" condition (see text) for Uprobe ~ 1 
meV. The intermediate dashed line in panel (a) represents the 
"dynamical instability" condition for uiprobe = 1 AteV. The ex- 
perimental points for T° (stars) and for T* measured with 
fast (crosses) and slow (squares) probes for LSCO are from 
Ref. those for Bi2212 are from Refs. The critical 

temperatures Tc are also shown (pluses). 

We stressed in several papers B that the actual onset 



of a fully developed CO phase is competing with local 
or coherent pair formation, which modifies the fermionic 
spectrum stabilizing the system against the electronic 
CO transition. There is the complementary possibil- 
ity that the charge-ordered state can never be reached 
because g is below gmin- The corresponding observable 
quantity identified by the experiments in Ref. pl[ | is the 
critical microstrain Sc, below which no stripe phase can 
be observed. Therefore the cuprates can deviate from 
CO criticality not only hy x — Xc in the (overdoped) 
quantum-disordered region, and by T or T — Tco{x) 
in the quantum-critical or underdoped region, but also 
because the microstrain e is smaller or larger than Ec, 
thereby tending to a homogeneous or an inhomogeneous 
phase respectively. As long as e — ec is not too large and 
T is finite, for x Xc the quantum-critical region, with 
strong dynamical fiuctuations, can be reached in both 
cases. According to the proposed mechanism of pairing 
mediated by critical fluctuations, the stronger are these 
fluctuations and the larger is T^. Indeed in the experi- 
ment the mercury compound IIgl212, having e ^ £c, has 
the the largest while IIgl201 (with e < Ec) and Bi2212 
or LSCO (with e > Ec) have a lower Tc- Moreover, since 
e < in IIgl201, we expect no Tco{x), and therefore 
no r*, for this material. 



B. Charge-ordering transition at finite temperature 
and pseudogap formation 

In this section we calculate the temperature for the 
onset of the CO instability within the mean-field approx- 
imation, T^q{x), and its corrected value, Tco{^)j in the 
presence of the leading fluctuations beyond mean field. 
For establishing the onset of the CO phase coming from 
high-doping and high-temperature regime, the direct in- 
volvement of the spin degrees of freedom can be safely 
ignored. 

For T < T^Q the CO fluctuations become substan- 
tial, leading to a reduction of the quasiparticle density of 
states. We therefore identify our mean-fleld line T^q{x) 
with the weak-pseudogap crossover line T''\x) observed 
in Knight-shift, transport, and static susceptibility mea- 
surements |l) , as the incipient reduction of the quasipar- 
ticle density of states. The fluctuation-corrected critical 
line Tco{x) T^q{x)], where the correlation length ^ 
diverges, and the scattering (^ becomes effectively singu- 
lar, is instead closely tracked by the pseudogap crossover 
line T*{x), according to the discussion of Sec. II. 

All the complicated formal structure of the quasi- 
particle scattering, mediated by phonons and by 
the Coulomb interaction within the slave-boson ap- 
proach, is, at the end, simply represented, at least 
near criticality, by a RPA resummation r(q, ciJn) w 
T4//(q)/[l + K//(q)n(q, w„)] of an effective static in- 
teraction Veff{(\) = f/(q) -I- Vb(q) - At pi. Here H is 
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the fermionic polarization bubble, A = 2g'^/tu)Q is the 
dimensionless e-ph coupling, and ?7(q) ~ A + B\(\\^ is 
the residual short-range repulsion between quasiparticles, 
and Vc'(q) is the Fourier transform of the e-e Coulomb 
repulsion. For the correspondence of Ve//(q) with the 
parameters of the original HH model see Ref. . The 
mean-field CO instability condition, which occurs for rea- 
sonable values A ~ 1, is 1 -|- T4//(qc)n(qc, w = 0) = 0. 
At T = this determines qc and the position of the 
mean- field QCP x^. By expanding 1 + Vez/H near the 



instability at T = 0, we find 



oc X 



At T > 0, 



considering the T dependence of the bare polarization 
bubble n, the instability condition determines the mean- 
field critical line T^q{x) (the solid- line curves in Fig. 
2), which ends at a;°. We evaluate Tqq{x) by taking 
standard quasiparticle (i.e., dressed by the slave bosons) 
band parameters {t — 200 meV, = 70 meV, leading 
to A = 200 meV and i? = 170 meV in the residual repul- 
sion U). Vc — 220 meV and g = 210 meV are adjusted 
to match with the extrapolated experimental values of 
T^{x) for LSCO § [x^ « 0.22) with the mean-field QCP 
x^. Without any further adjustment of the parameters 
the curve T^q {x) agrees remarkably well with the exper- 
imental data for T^(x). Similar parameters are taken for 
Bi2212 to fit the data in Refs. (^[|2|], with Vc ^ 220 eV 
and g = 230 eV. 




FIG. 3. The vertex and selfenergy corrections beyond RPA 
for the fermionic bubbles (solid line) due to critical CO fluc- 
tuations [the wavy line represents the singular part of the 
effective interaction, Eq. (nl)]. 



coupling constant u ~ V/t^. Since the phonons are re- 
sponsible for the instability, loq appears as the frequency 
cut-off jlll. The fluctuation-corrected instability condi- 
tion is = 0. As it stands, Eq. (^ is written for the 
two-dimensional Cu-0 planes. At T = 0, Eq. (^) with 
= leads to a finite shift of the two-dimensional (2D) 
QCP x^- Xc oc u)o (see Fig. 2). At finite T, purely 2D 
fluctuations suppress the transition. However, by con- 
sidering the more realistic anisotropic 3D character of 
the critical fluctuations, we obtain a finite fluctuation- 
corrected transition line Tcoi^): which is reported as 
the lowest dashed line in Fig. 2(a). Thus, the inclusion 
of fluctuations brings the critical line from temperatures 
of the order of typical electronic energies {T^q ^ down 
to much lower temperatures Tco of the order of the ob- 
served T*'s. 

The spread in the measured values of T* depending on 
the typical frequency of the experimental probe ujprobe 
(see, e.g., Ref. j|l|) suggests, however, within our identi- 
fication of Tco with T*, that the CO instability may be 
"dynamical". In this case, the critical condition = 
is replaced by = 1^ probe in the self-consistency con- 
dition, Eq. (||). Two examples are reported in Fig. 2 for 
t^probe — 1 meV (dot-dashed fine in Figs. 2(a,b), typical 
of neutron scattering experiments, and Wprobe = 1 A'eV 
(second dashed line from bottom in Fig. 2(a), typical 
of static experiments (NQR, NMR). The corresponding 
experimental data for T* in LSCO are also reported for 
comparison. The coupling V between the charge fluctu- 
ations and the quasiparticles, which is the only parame- 
ter for which an a priori estimate is difficult, is fixed by 
matching the position of the fluctuation-corrected QCP 
with the T — extrapolation of the experimental val- 
ues of T*. We used 7 = 0.7 eV"! and 7 = 0.4 eV'^ 
for LSCO and Bi2212 respectively, and V = 0.54 eV for 
both. Note that, similarly to the case of Tqq{x), the 
agreement between the calculated Tcoi^) and the ex- 
perimental points for T* is obtained without adjustment 
of other parameters. 



Due to the strong anisotropy of the layered cuprates, 
the contribution of the fluctuations related to any insta- 
bility is significant |^,^. Specifically the CO fluctua- 
tions which mediate the critical effective interaction, Eq. 
(|l|), can be included in the polarization bubble If (en- 
tering the instability condition) via the diagrams of Fig. 
3, leading to corrections beyond RPA. From the explicit 
evaluation of these diagrams we flnd the self-consistent 
correction to the inverse square correlation length 

V 



UuT 



E 



^-2 + |q-q,|2+7|c^„|^ 



(2) 



where a;„ are the bosonic Matsubara frequencies, and the 



III. NOVEL ISOTOPE EFFECTS 

In the present framework, any mechanism shifting the 
position of the QCP and the line Tcoix) is mirrored 
by a corresponding shift of the superconducting critical 
line Tc{x) and of the pseudogap formation temperature 
T*{x). This is the case of the shift induced by isotopic 
substitution, and we want to connect our results with 
the observation of isotopic effects on Tc [ p8|j29| and on 
T* [Eol-^. As shown in Sec. II A, the mean-fleld weak- 



pscudogap crossover temperature Tqq ^ depends on 
the phonon frequency only through the dimensionless e- 
ph coupling A, which is known to display no isotopic de- 
pendence. Therefore T^q, and consequently and the 
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portion oiTc{x) near and above T^{x), are not expected 
to display isotope effects. 

On the other hand, the CO fluctuations, as shown 
above, crucially involve loq, and new physical effects can 
be induced by the isotopic substitutions, on fluctuation- 
dependent physical quantities, such as Xc and Tcoix). 
Since uq decreases with increasing ionic mass, the effect 
of fluctuations is reduced, and Xc and Tco{x) remain 
closer to their mean-field values. In Fig. 4, we report 
the line Tcoix), calculated via Eq. (||), with the same 
parameters used in Sec. II B to fit the T* data of LSCO, 
together with its isotopic shift calculated for ^^O ^^O 
substitution (i.e., for a five percent reduction of ujq). 
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FIG. 4. Calculated effect of tlie isotopic cliange ^''O O 
(i.e. LOO = TO meV and loo = 66 meV) on T* in LSCO, botli 
for a fast- and a slow-probe measurement. Tfie inferred shift 
on Tc is also reported, although hardly visible on this scale. 

Contrary to standard theories based on CO pseudo- 
gap [Q, the isotopic shift is positive (i.e. anomalous) 
for Tco ~ T*. Moreover, if the slope of Tco{x) is 
large, a small isotopic shift in Xc can result in a sub- 
stantial shift in Tco ~ T*, i.e., the steeper is T* , the 
larger is the isotopic effect. Since the slope of Tco{x) 
increases by increasing ujprobe, whereas the QCP is un- 
shifted, the isotope effect on Tco ^ T* is enhanced, and 
we have the general trend that faster probes should de- 
tect a larger isotope effect on T*. Although we cannot 
account for the near-absent or negative isotope effect on 
T* within the almost static probes in YBa2Cu408 pO| , ^ , 
this general trend is in qualitative agreement with a much 
stronger effect observed in the isostructural HoBa2Cu408 
(HBCO-124) with fast neutron scattering Indeed, 
this fast-probe experiment should be represented by the 
curve Tco{x) corresponding to LOprobe = 1 nieV in Fig. 
4. The huge isotope effect on T* observed in IIBCO-124 
suggests that the curve Tco{x) in this system is steeper 
than in LSCO. 

As far as in the underdoped region is concerned, 
we can only infer some expected consequences. Assum- 
ing that the shift in Xc produces a rigid shift of Tc{x) 
along the doping axis, the isotopic effect for Tc is nega- 



tive (i.e. normal) upon reducing wq, and much smaller 
than for T* (the Tc{x) curve being much flatter), in 
agreement with long-standing experiments ||2^,^. This 
large difference in the isotope effect for Tc and T*, 
(ATc/AM) / (AT* /AM) < 1, is indeed experimentally 
observed in IIBCO-124, and reported in ReL |Q, where 
a "striking similarity between isotopic substitution and 
underdoping with respect to both Tc and T*" is pointed 
out. Although we are not aware of any systematic analy- 
sis of the doping dependencies of Tc, T* and their isotopic 
shifts in IIBCO-124, allowing for a strict comparison, 
this observation finds its natural interpretation within 
our theory, where the isotopic substitution increases Xc 
and is therefore nearly equivalent to underdoping. 

In the strongly underdoped materials, where substan- 
tial isotope effects in the penetration depth and in 
the x-ray absorption have been observed, our results 
can provide only qualitative indications, as other effects 
(magnetic, polaronic, lattice-pinning) not included in our 
HH model become relevant. 

IV. CONCLUSIONS 

In conclusion we have revisited the CO scenario, and 
explicitly evaluated the mean-field and the fluctuation- 
corrected critical lines, T^q(x) and Tco{x), for the onset 
of the charge inhomogeneous phase, which should rep- 
resent the stripe phase of the cuprates, when present. 
The nearly singular effective interaction among the elec- 
trons mediated by the CO fluctuations near the insta- 
bility, suggests the natural identification of the experi- 
mentally observed weak-pseudogap crossover line T°(x) 
with T^q(x), and of the pseudogap formation temper- 
ature T* with Tcoix)- We obtain in this way a good 
fitting of the experimental data, and by assuming that 
the CO transition has a dynamical character, we account 
for the spread of the measurements of T*(x), depend- 
ing on the characteristic time scale of the experimental 
probe. The existence of the CO transition, depending 
on the strength of the e-ph coupling, as produced by our 
zero-temperature phase diagram, can give an answer to 
the puzzling question why features of formed stripes are 
present in some cuprates and absent in others. The CO 
instability condition at mean-field level depends only on 
the dimensionless e-ph coupling which is not changed by 
isotopic substitution. Therefore, in our model, T^pix) 
should not undergo any isotopic shift. On the contrary, 
Tcoix), and therefore T*, arc strongly modified by the 
effect of CO fluctuations, which in our approach is pro- 
portional to the phonon frequency loq. Within our model, 
therefore, T* shows a novel isotope effect produced by 
charge fluctuations, with a positive shift, which depends 
on the slope of Tcoix)- It results, therefore, that the iso- 
tope effect is stronger when detected with faster probes. 
We also infer that Tc(a;) should have no isotope effect in 
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the optimal and overdoped regimes, and a small negative 
effect in the underdoped region, which increases upon un- 
derdoping. In this new isotope effect, the phonon do not 
appear directly as mediators of pairing, but are indirectly 
involved, via the CO fluctuations, and qualitatively ex- 
plain the complex behavior of T* and upon isotopic 
substitution. 
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